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I review arguments for the existence of a critical point E in the QCD phase diagram 
as a function of temperature T and baryon chemical potential fi. I describe how heavy 
ion collision experiments at the SPS and RHIC can discover the tell-tale signatures of 
such a critical point, thus mapping this region of the QCD phase diagram. I contrast 
the different ways in which the matter produced in a heavy ion collision can be driven 
out of equilibrium: quenching out of equilibrium (possible, but not guaranteed, if the 
transition region is traversed at /i <C he) vs. slowing out of equilibrium (guaranteed for 
M ~ He) vs. bubbling out of equilibrium (possible, but not guaranteed, for ji ^> fi E ). 
Quenching or bubbling create and amplify distinct, detectable, non-gaussian fluctuations. 
In contrast, slowing out of equilibrium reduces the magnitude of the specific, detectable, 
gaussian fluctuations which signal the presence of the critical point. 

1. The Critical Point 

One goal of relativistic heavy ion collision experiments is to explore and map the QCD 
phase diagram as a function of temperature and baryon chemical potential. Recent theo- 
retical developments suggest that a key qualitative feature, namely a critical point which 
in a sense defines the landscape to be mapped, may be within reach of discovery and 
analysis by the CERN SPS or by RHIC, as data is taken at several different energies 
. The discovery of the critical point would in a stroke transform the map of the QCD 
phase diagram from one based only on reasonable inference from universality, lattice gauge 
theory and models into one with a solid experimental basis Qj. 

In QCD with two massless quarks (m u ^ = 0; m s = oo) the phase transition at which 
chiral symmetry is restored as T is increased with fi = is likely second order and 
belongs to the universality class of 0(4) spin models in three dimensions Q. Below T c , 
chiral symmetry is broken and there are three massless pions. At T = T c , there are four 
massless degrees of freedom: the pions and the sigma. Above T = T c , the pion and sigma 
correlation lengths are degenerate and finite. 
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Figure 1. Sketch of the QCD phase diagram as a function of temperature T and baryon 
chemical potential \i. Chiral symmetry is broken at low T and fi. As T is increased, chiral 
symmetry is approximately restored via a smooth crossover to the left of E or a first order 
phase transition to the right of E. The symmetry is only approximately restored because 
the light quarks are not massless. At the critical point E at which the line of first order 
phase transitions ends, the transition is second order and is in the Ising universality class. 
(At large [i and small T, there are color superconducting phases.) The Ising model r-axis 
and h-axis and the trajectories a, b and c will be discussed below. 



In nature, the light quarks are not massless. Because of this explicit chiral symmetry 
breaking, the second order phase transition is replaced by an analytical crossover: physics 
changes dramatically but smoothly in the crossover region, and no correlation length 
diverges. This picture is consistent with present lattice simulations ||, which suggest 
T c ~ 140 - 190 MeV [§. 

Arguments based on a variety of models [ff|-[T4]1 indicate that the transition as a function 
of T is first order at large This suggests that the phase diagram features a critical 
point E at which the line of first order phase transitions present for fi > fiE ends, as 
shown in Figure 1.0 At fiE, the phase transition at T = Te is second order and is in 
the Ising universality class [|11|,|T^]. Although the pions remain massive, the correlation 



length in the a channel diverges due to universal long wavelength fluctuations of the order 
parameter. This results in characteristic signatures, analogues of critical opalescence in 
the sense that they are unique to collisions which freeze out near the critical point, which 
can be used to discover E [JI],^. 

The position of the critical point is, of course, not universal. Furthermore, it is sensitive 
to the value of the strange quark mass, he decreases as m s is decreased Q, and at some 
m c s , it reaches he — and the transition becomes entirely first order |T5|]. The value of 
m c s is an open question, but lattice simulations suggest that it is about half the physical 
strange quark mass [|H],[T7]], although these results are not yet conclusive [[ 



Of course, 



2 The realization that this phase transition is not simply a transition at which the chiral condensate 
vanishes, but is instead a transition at which a phase with nonzero chiral condensate and a phase with 
nonzero superconducting condensate compete and coexist has further strengthened the arguments in 
favor of a first order phase transition on the /x-axis in two-flavor QCD ||[ll 13 14j ]. 
If the up and down quarks were massless, E would be a tricritical point, at which the first order transition 
becomes second order. 
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experimentalists cannot vary m s . They can, however, vary fi. The AGS, with beam 
energy 11 AGeV corresponding to y/s = 5 GeV, creates fireballs which freeze out near 
H ~ 500 - 600 MeV @. When the SPS runs with ^ = 17 GeV (beam energy 158 
AGeV), it creates fireballs which freeze out near ii ~ 200 MeV ]nj. RHIC will make 
even smaller values of \x accessible. By dialing y/s and thus /a, experimenters can find the 
critical point E. Finding E would confirm that m c s is less than the physical value of m s , 
and that the phase transition at lower chemical potentials (// < is a smooth crossover. 

2. Discovering the Critical Point 

Predicting fi E , and thus suggesting the y/s to use to find E, is beyond the reach of 
present theoretical methods because \ie is both nonuniversal and sensitively dependent 
on the mass of the strange quark. Crude models suggest that lie could be ~ 600 — 800 
MeV in the absence of the strange quark |TT],|12j; this in turn suggests that in nature /i£ 



may have of order half this value, and may therefore be accessible at the SPS if the SPS 
runs with y/s < 17 GeV. However, at present theorists cannot predict the value of \xe 
even to within a factor of two. The SPS can search a significant fraction of the parameter 
space; if it does not find E, it will then be up to the RHIC experiments to map the 
He < 200 MeV region. 

Locating E on the phase diagram can only be done convincingly by an experimental 
discovery. Theorists can, however, do reasonably well at describing the phenomena that 
occur near E, thus enabling experimenters to locate it. This is the goal of Ref. [0. 
The signatures proposed there are based on the fact that E is a genuine thermodynamic 
singularity at which susceptibilities diverge and the order parameter fluctuates on long 
wavelengths. The resulting signatures are nonmonotonic as a function of y/s: as this 
control parameter is varied, we should see the signatures strengthen and then weaken 
again as the critical point is approached and then passed. 

The simplest observables to use are the event-by-event fluctuations of the mean trans- 
verse momentum of the charged particles in an event, pt, and of the total charged multi- 
plicity in an event, N. The fluctuations observed in SPS collisions by NA49 are as perfect 
gaussians as the data statistics allow, as expected for freeze-out from a system in thermal 
equilibrium [pOj] . The data on multiplicity fluctuations show evidence for a nonthermo- 
dynamic contribution, which is to be expected since the extensive quantity N is sensitive 
to the initial size of the system and thus to nonthermodynamic effects like variation in 
impact parameter. The contribution of such effects to the fluctuations have now been es- 
timated pl,p2|; the combined thermodynamic and nonthermodynamic fluctuations are in 



satisfactory agreement with the data ||22|| . The width of the event-by-event distribution of 
the mean pt of the charged pions in a single event[] is in good agreement with predictions 
based on noncritical thermodynamic fluctuations @. That is, NA49 data are consistent 
with the hypothesis that almost all the observed event-by-event fluctuation in mean p?, 
an intensive quantity, is thermodynamic in origin. This bodes well for the detectability 
of systematic changes in thermodynamic fluctuations near E. 



4 This width can be measured even if one observes only two pions per event |23| ]; large acceptance data as 
from NA49 is required in order to learn that the distribution is gaussian, that thermodynamic predictions 
may be valid, and that the width is therefore the only interesting quantity to measure. 
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One analysis described in detail in Ref. @ is based on the ratio of the width of the true 
event-by-event distribution of the mean to the width of the distribution in a sample 
of mixed events. This ratio was called y/~F. NA49 has measured y/F = 1.002 ± 0.002 
[[20|,^j , which is consistent with expectations for noncritical thermodynamic fluctuations.^ 
Critical fluctuations of the a field, i.e. the characteristic long wavelength fluctuations 
of the order parameter near E, influence pion momenta via the (large) ann coupling 
and increase \/F @. The effect is proportional to £f ree zeout> where £f re ezeout is the cr-field 
correlation length of the long- wavelength fluctuations at freezeout 0. If £f ree zeout ~ 3 fm, 
the ratio y~F increases by ~ 3 — 5%, ten to twenty times the statistical error in the present 
measurement ||. This observable is valuable because data on it has been analyzed and 
presented by NA49, and it can therefore be used to learn that Pb+Pb collisions at 158 
AGeV do not freeze out near E. 

Once E is located, however, other observables which are more sensitive to critical effects 
will be more useful. For example, a \AF so ft, defined using only the softest 10% of the pions 
in each event, will be much more sensitive to the critical long wavelength fluctuations. The 
higher pr pions are less affected by the a fluctuations [Q] , and these relatively unaffected 
pions dominate the mean pt of all the pions in the event. This is why the increase in \J~F~ 
near the critical point will be much less than that of \J F so f t . 

The multiplicity of soft pions is an example of an observable which may be used to 
detect the critical fluctuations without an event-by-event analysis. The post-freezeout 
decay of sigma mesons, which are copious and light at freezeout near E and which decay 
subsequently when their mass increases above twice the pion mass, should result in a 
population of pions with pt ~ m n /2 which appears only for freezeout near the critical 
point 0. If ^freezeout > this population of unusually low momentum pions will 

be comparable in number to that of the "direct" pions (i.e. those which were pions at 
freezeout) and will result in a large signature. 

The variety of observables which should all vary nonmonotonically with y/s (and should 
all peak at the same y/s) is sufficiently great that if it were to turn out that //e < 200 
MeV, making E inaccessible to the SPS, all four RHIC experiments could play a role in 
the study of the critical point. 



3. Slowing Out of Equilibrium: How Large Can £ Grow? 



The purpose of Ref. is to estimate how large £f ree zeout can become, thus making the 
predictions of Ref. @ for the magnitude of various signatures more quantitative. In an 
ideal system of infinite size which was held at T = T#; /j = /Zg for an infinite time, the 
correlation length £ would be infinite. Ref. M estimated that finite size effects limit £ 



to be about 6 fm at most. We argue in Ref. []53| that limitations imposed by the finite 
duration of a heavy ion collision are more severe, preventing £ from growing larger than 
about 2/T E ~ 3 fm. 



5 In an infinite system made of classical particles which is in thermal equilibrium, ^/F = 1. Bose effects 
increase y/F by 1 — 2% [p4|]2]]; an anticorrelation introduced by energy conservation in a finite system 
— when one mode fluctuates up it is more likely for other modes to fluctuate down — decreases \fF by 
1 — 2% j2|; two-track resolution also decreases y/F by 1 — 2% ^0|. The contributions due to correlations 
introduced by resonance decays and due to fluctuations in the flow velocity are each much smaller than 

i%0. 
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The nonequilibrium dynamics analyzed in Ref. [25]| is guaranteed to occur in a heavy 



ion collision which passes near E, even if local thermal equilibrium is achieved at a higher 
temperature during the earlier evolution of the plasma created in the collision. We assume 
early thermal (although not necessarily chemical) equilibration, and ask how the system 
evolves out of equilibrium as it passes E. For the present, assume that the system cools 
through the critical point E as sketched in trajectory (b) of Figure 1. If it were to cool 
arbitrarily slowly, £(T) = £ eq (^) would be maintained at all temperatures, and £ would 
diverge at T%. However, it would take an infinite time for £ to grow infinitely large. Indeed, 
near a critical point, the longer the correlation length, the longer the equilibration time, 
and the slower the correlation length can grow. This critical slowing down means that the 
correlation length cannot grow as fast as £ eq , and the system cannot stay in equilibrium. 

We use the theory of dynamical critical phenomena to describe the effects of critical 
slowing down on the time development of the correlation length £(£), using the three- 



dimensional Ising model to describe £ eq (T, //) near E [p5| . In the Ising model, the order 
parameter and the correlation length are functions of the reduced temperature r and the 
magnetic field h. (In the Ising model, r is defined as (T — T c )/T c and is usually called 
£; we reserve the symbol t for time, however.) We assume that the r- and /i-axes are 
as sketched in Figure 1, although in reality they will be deformed. We consider cooling 
along trajectories that pass through E at various angles, like (a) and (b) sketched in 
Figure 1 and also consider trajectories like (c) which come near to but miss E. Our 
results depend on the universal function describing £ e q(?", h), on the universal dynamical 
exponent z describing critical slowing down (perturbations away from equilibrium relax 
toward equilibrium on a timescale which scales with £ like A^ z ||26||), on the nonuniversal 
constant A, the nonuniversal constants which relate r and h to (T — T#) and (// — 
on Te which we take to be ~ 140 MeV, and finally on the cooling rate \dT/dt\ which we 
estimate to be 4 MeV/fm [g7|-|2jj[2§ . 



Our results have a number of consequences, detailed in the remainder of this section, 
which should be taken into account both in planning experimental searches for the QCD 
critical point, and in planning future theoretical work. 

We do indeed find that because of the critical slowing down of the long wavelength 
dynamics near E, the correlation length does not have time to grow as large as it would 
in equilibrium: we find £f re ezeout ~ 2/Tg ~ 3 fm for trajectories passing near E. 

Although critical slowing down hinders the growth of £, it also slows the decrease of £ 
as the system continues to cool below the critical point. As a result, £ does not decrease 
significantly between the phase transition and freezeout. 

Our estimate that £ does not grow larger than 2/Te is robust in three senses. First, it 
depends very little on the angle with which the trajectory passes through E. Second, it 
turns out to depend on only one combination of all the nonuniversal quantities which play 
a role. We call this parameter a; it is proportional to \dT/dt\~ l . Third, our results do not 

depend sensitively on a. We show that the maximum value of £ scales like a 1 + z "V< 3 < 5 m a 0,215 
p5[ .[] Thus, for example, \dT/dt\ would have to be a factor of 25 smaller than we estimate 



6 A scaling law of this form (but of course with different exponents) relating the maximum correlation 
length which is reached to the cooling rate with which a second order phase transition is traversed was first 
discovered in the theory of defect formation at a second order phase transition f3(i|| . In this context, the 
maximum correlation length reached during the phase transition sets the scale for the initial separation 
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in order for £ to grow to instead of 2/Tg. Although our results are robust in this 

sense, they cannot be treated as precise because our assumption that the dynamics of £ in 
QCD is described by the universal classical dynamics of the three-dimensional Ising model 
only becomes precise if £ ^> 1/Te, while our central result is that £ does not grow beyond 
~ 2/Te. A 3 + 1- dimensional quantum field theoretical treatment of the interplay between 
cooling and the dynamics of critical slowing down is not yet available, but promising first 
steps in this direction can be found in Ref. []34] . 



A result which is of great importance in the planning of experimental searches is that 
one need not hit E precisely in order to find it. Our analysis of trajectories like (c) 
of Figure 1 demonstrates that if one were to do a scan with collisions at many finely 
spaced values of the energy and thus /i, one would see signatures of E with approximately 
the same magnitude over a broad range of /x. The magnitude of the signatures will not 
be narrowly peaked as \i is varied. As long as one gets close enough to E that the 
equilibrium correlation length is (2 — 3)/Tg, the actual correlation length £ will grow to 
~ 2/Te- There is no advantage to getting closer to E, because critical slowing down 
prevents £ from getting much larger even if £ eq does. Data at many finely spaced values 
of \i is not called for.[] 

Knowing that we are looking for £f reeze out ~ 3 fm is very helpful in suggesting how to 
employ the signatures described in detail in Ref. 0. The excess of pions with p T ~ m^/2 
arising from post-freezeout decay of sigmas is large as long as £f ree zeout ~ l/ m 7n and 
does not increase much further if £f ree zeout is longer. This makes it an ideal signature. 
The increase in the event-by-event fluctuations in the mean transverse momentum of the 
charged pions in an event (described by the ratio \fF of Ref. @) is proportional to 
£freezeout- The results of Ref. suggest that for £f ree zeout ~ 3 fm, this will be a 3 — 5% 
effect. This is ten to twenty times larger than the statistical error in the present NA49 
data, but not so large as to make one confident of using this alone as a signature for E. 
The solution is to use signatures which focus on the event-by-event fluctuations of only 
the low momentum pions. Unusual event-by-event fluctuations in the pion momenta arise 
via the coupling between the pions and the sigma order parameter which, at freezeout, 
is fluctuating with correlation length £f ree zeout- This interaction has the largest effect on 
the softest pions 0. V-^soft, described in the previous section, is a good example of an 
observable which takes advantage of this. Depending on the details of the cuts used to 
define it, it should be enhanced by many tens of percent in collisions passing near E. Ref. 
[@] suggests other such observables, and more can surely be found. Together, the excess 
multiplicity at low momentum (due to post-freezeout sigma decays) and the excess event- 
by-event fluctuation of the momenta of the low momentum pions (due to their coupling to 



between defects — vortices in a superfiuid or in liquid crystals, for example — created as the system cools 
through the transition. This initial network of defects coarsens at later times. The Ising phase transition 
of interest to us creates no defects. It is nevertheless very pleasing that scaling laws analogous to the one 
we need have been tested quantitatively in numerical simulations of defect formation and dynamics |31[ ] 
and, furthermore, are supported by data from experiments on liquid crystals p2| and superfiuid 3 He ]35( |. 
7 Analysis within the toy model of Ref. suggests that in the absence of the strange quark, the range 
of fi over which £ oq > 2 fm is about A/i ~ 120 MeV for he ~ 800 MeV. Similar results can be obtained 
within a random matrix model |]l2|. It is likely over-optimistic to estimate A/i ~ 120 MeV when the 
effects of the strange quark are included and fiE itself is reduced. A conservative estimate would be to 
use the models to estimate that Afi/fiE ~ 15%. 
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the order parameter which is fluctuating with correlation length £f ree zeout) should allow a 
convincing detection of the critical point E. Both should behave nonmonotonically as the 
collision energy, and hence //, are varied. Both should peak for those heavy ion collisions 
which freeze out near E, with £f ree zcout ~ 3 fm. 



4. Quenching Out of Equilibrium vs. Slowing Out of Equilibrium vs. Bubbling 
Out of Equilibrium 

As the matter created in a heavy ion collision cools through the QCD phase transition, 
it may be driven out of equilibrium. The way in which this may happen is quite different 
depending on whether the phase transition is traversed with \i 3> \i Ey /i C /i E or /i ~ he- 

If /i > the phase transition is first order; if (/i — n E ) is large enough, the transition 
may be sufficiently strongly first order that it proceeds via the nucleation of well-separated 
bubbles of hadronic matter, which then grow. This bubbling results in a spatially inho- 
mogeneous hadronic system. If the system freezes out before it gets re-homogenized, the 
result will be large, non-gaussian fluctuations in hadron multiplicity as a function of rapid- 
ity [Q. These overdense and underdense regions in rapidity would provide a distinctive 
signature, and have not been seen. 

If /i <C fiE, the phase transition is a smooth crossover. This is the region of the diagram 
where the transition may be traversed most quickly, with the most rapid cooling rate. 
(For /i ~ he, the specific heat is unusually large and therefore even though the rate of 
decrease of energy density is not unusual, the cooling rate \dT/dt\ is unusually small. 
For ^> fiE, neither of the two phases which coexist at the first order phase transition 
has an unusually large specific heat. The system nevertheless spends a long time at the 
coexistence temperature T c because of the release of latent heat. Thus, \dT/dt\ is allowed 
to be largest for fi <C fiE-) Although the smooth crossover with nondivergent specific heat 
allows a large \dT/dt\ in principle, the cooling rate which is actually achieved will depend 
on the system size and on the duration of time between the initial collision and the phase 
transition. If conditions can be found in which \dT/dt\ can be made large enough, the 
dynamics may then be well-modelled as a quench, in which most of the degrees of freedom 
in the system are imagined to cool arbitrarily rapidly. In this circumstance, the long 
wavelength modes of the pion field become unstable to exponential growth pTj, leading 



to the formation of long wavelength pion oscillations. Once excited, these oscillations may 
be further amplified by parametric resonance These long wavelength disorientations 



of the chiral condensate result in characteristic fluctuations in the ratio of the number of 
low pt neutral pions to the number of low px charged pions |3T||3"T| , |4Ti|] . Other signatures 



of DCC oscillations have also been discussed gTJ. WA98 [§g] and NA49 [£(J have looked 
for signatures of DCC oscillations in 158 AGeV PbPb collisions, but none have been seen. 

The considerations above suggest that RHIC collisions provide more favorable condi- 
tions for DCC production than SPS collisions, because /i is smaller and one is more likely 
to have a smooth crossover which allows large \dT/dt\ . However, the hadronic systems 
produced in RHIC collisions will be large, will traverse the phase transition late, and 
will freeze out even later. The late transition reduces the cooling rate, making DCC 



8 As seems indicated 38 by the long period of time over which amplification is seen in the simulations of 
Ref. |7). 
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formation less likely, and the late freeze out makes wash-out of DCC signatures a worry. 
This suggests that the most favorable conditions in which to look for DCC signatures 
may be RHIC collisions of medium-sized ions, or AuAu collisions at RHIC with impact 
parameters of order the nuclear radius. (That is, midway between central and peripheral 
collisions.) Quite different considerations have lead Asakawa, Minakata and Miiller to 
propose that non-central RHIC collisions provide the most favorable conditions for DCC 
formation fll3[j : these collisions feature strong electromagnetic E • B fields which may, 
through the anomaly, result in a coherent excitation of the 7r° field. There are therefore 
three reasons to suggest that non-central RHIC collisions may yield DCC signatures: the 
smooth crossover and small system size allow the system to cool through the transition 
faster, favoring the instability driven growth of DCCs; the electromagnetic kick may fur- 
ther drive their growth; and, the smaller system size will hasten freezeout. The first and 
third reasons also apply to RHIC collisions of medium-sized ions. 

For /i ~ fiE, the specific heat is larger in the vicinity of the critical point than any- 
where else on the phase diagram, and the cooling rate \dT/dt\ is correspondingly slow. 
This slow cooling means that those collisions which create matter which cools through 
the transition region with [i ~ [i^ are least likely to be driven far from equilibrium. That 
is, we do not expect DCC signatures or signatures of bubbling in these collisions. The 
interesting twist described in the previous section is that it is nevertheless true that such 
collisions cannot stay precisely in equilibrium [pj]. Even though the system cools through 
the vicinity of E slowly, the correlation length still grows much too slowly to become as 
long as it would in equilibrium. We have discussed this phenomenon of slowing out of 
equilibrium at length, as it determines the magnitude of signatures of the critical point. 
Here, we note that this non-equilibrium physics is much milder than quenching or bub- 
bling. Quenching or bubbling create more fluctuations than would be obtained in thermal 
equilibrium. They also create fluctuations which are non-gaussian in characteristic ways. 
In contrast, slowing out of equilibrium serves to reduce the growth of fluctuations relative 
to what would be achieved if equilibrium were maintained. The critical point E will be 
signalled by an increase in the magnitude of suitably chosen event-by-event fluctuations, 
for example parametrized by the growth in yJ~F and y/F so{t , but these fluctuations are 
expected to remain gaussian. Slowing out of equilibrium means that these gaussian fluc- 
tuations are smaller than they would be if equilibrium were maintained near E. They are 
nevertheless larger than the equilibrium thermal fluctuations far from E and hence result 
in nonmonotonic signatures which peak for collisions which freezeout near E^\ 

We have learned much from the beautiful gaussian event-by-event fluctuations observed 
by NA49. The magnitude of these flucutations are consistent with the hypothesis that 
the hadronic system at freezeout is in approximate thermal equilibrium. These and other 
data show no signs that the system at freezeout has recently been bubbled or quenched. 
There is also no sign of the larger, but still gaussian, fluctuations which would signal 
freezeout near the critical point E. Combining these observations with the observation 



9 In addition to restricting the growth of fluctuations, critical slowing down retards the decrease of the 
fluctuations after the transition. In addition, the large specific heat and consequent slow cooling rate 
result in an elevated freezeout temperature Q|. This suggests that signatures of E will not be washed 
out before freezeout; if necessary, however, smaller ions or non-central collisions can be used to further 
elevate the freezeout temperature fi^H. 
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of several tantalizing indications that the matter created in SPS collisions is not well 
described at early times by hadronic models [|5| suggests that collisions at the SPS may 
be exploring the crossover region to the left of the critical point E, in which the matter 
is not well-described as a hadron gas but is also not well-described as a quark-gluon 
plasma. This speculation could be confirmed in two ways. First, if the SPS is probing 
the crossover region then the coming experiments at RHIC may discover direct signatures 
of an early partonic phase, which are well-described by theoretical calculations beginning 
from an equilibrated quark-gluon plasma. Second, if 158 AGeV collisions are probing 
the crossover region not far to the left of the critical point E, then running at lower 
energies would result in the discovery of E. If, instead, RHIC were to discover E with 
He < 200 MeV, that would indicate that the SPS experiments have probed the weakly 
first order region just to the right of E. Regardless, discovering E would take all the 
speculation out of mapping this part of the QCD phase diagram. 
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